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ABSTRACT 

The Total Irradiance Monitor (TIM) from NASA’s SOlar Radiation and Climate Experiment can detect changes in 
the total solar irradiance (TSI) to a precision of 2 ppm, allowing observations of variations due to the largest X-class 
solar flares for the first time. Presented here is a robust algorithm for determining the radiative output in the TIM 
TSI measurements, in both the impulsive and gradual phases, for the four solar flares presented in Woods et al., as 
well as an additional flare measured on 2006 December 6. The radiative outputs for both phases of these five flares 
are then compared to the vacuum ultraviolet (VUV) irradiance output from the Flare Irradiance Spectral Model 
(FISM) in order to derive an empirical relationship between the FISM VUV model and the TIM TSI data output 
to estimate the TSI radiative output for eight other X-class flares. This model provides the basis for the bolometric 
energy estimates for the solar flares analyzed in the Emslie et al. study. 
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1. INTRODUCTION 

The total solar irradiance (TSI) is the energy released by 
the Sun across the entire electromagnetic spectrum. TSI con- 
tributions due to solar flares are normally masked by other 
solar variations, such as acoustic oscillations and granulation, 
which occur on timescales similar to the initial rapid bolometric 
enhancement from solar flares. In the last decade, large solar 
flares have been detected from the TSI signal directly from the 
Total Irradiance Monitor (TIM; Kopp et al. 2005b; Kopp & 
Lean 2011) and through an epoch superposition technique us- 
ing the Variability of solar IRradiance and Gravity Oscillations 
(VIRGO) experiment (Kretzschmar et al. 2010). Solar flares 
that exhibit measurable visible light enhancements have been 
deemed “white light flares” (Hudson et al. 1992). It is possible 
that these white light emissions contribute large amounts of en- 
ergy to the overall flare energy budget. Their contribution to the 
TSI energy budget is still an ongoing investigation. This paper 
aims to help constrain the variation in the TSI due to flares. 

Solar flares are the dominant cause of short timescale vacuum 
ultraviolet (VUV; 0.1-200 nm) emission. VUV wavelengths 
normally contribute a small portion to the TSI, but can vary by 
a large factor during solar activity. These emissions are one of 
the major drivers and energy input source of the Earth’s upper 
atmosphere. The increase in VUV irradiance and charged parti- 
cles from the Sun can affect Earth’s thermosphere, mesosphere, 
and ionosphere (Qian et al. 2011). Radio transmissions can be 
interrupted and global positioning systems can be disturbed. 
An increase in knowledge of the flare VUV and TSI energy 
release is desired for studies of comparing various energies in 
solar eruptive events (SEEs), such as was done in Emslie et al. 
( 2012 ). 

Physically, a solar flare is defined as an event where stored 
magnetic energy is released via an instantaneous (less than a 
second) rearrangement of stressed magnetic field lines, presum- 
ably magnetic reconnection, that accelerates charged particles 
both toward the solar surface and outward into space (Hudson 
2011). The radiative enhancement stems from an increase in 


particle flux in the solar atmosphere and can be partitioned into 
two phases, an impulsive phase (IP) and a gradual phase (GP). 
The initial flare energy release accelerates charged particles in- 
ward toward the more dense plasma by the solar surface. This 
energy input to the lower solar atmosphere heats the plasma and 
drives an increase of the hard X-ray flux, via bremsstrahlung. 
The heated particles confined within the magnetic loops subse- 
quently radiate thermally in soft X-rays and in some portions of 
the VUV. 

The initial hard X-ray temporal enhancement is called the 
IP and the GP is the radiative cooling of the heated plasma in 
soft X-rays and VUV. The IP generally lasts for 5-10 minutes, 
and the GP peaks shortly after the IP and can last for several 
hours. This observed relation between the hard and soft X-ray 
components, deemed the Neupert effect (Dennis & Zarro 1993; 
Veronig et al. 2002; Benz 2008; Hudson 201 1), originates from 
the work of Neupert (1968), who noticed that the integral of 
the centimeter radio flux roughly matched the soft X-ray flux. 
The acceleration of electrons outward to relativistic speeds 
generating radio emission is linked to the hard X-ray emission 
from the inward-penetrating electrons. Investigations of white 
light flares in spectral lines and continuum bands with multiple 
satellites demonstrate their temporal and spatial coherence with 
hard X-ray emission (Martinez Oliveros et al. 2011; Watanabe 
et al. 2013). This supports the case of white light emission 
occurring primarily during the IP. 

This study is focused on the radiative aspect of large X-class 
solar flares. An X-class flare is a flare that has a peak irradiance 
greater than 1 0 4 W m 2 in the 0. 1-0.8 nm channel of the 
Geostationary Operational Environmental Satellite (GOES) 
X-Ray Sensor (XRS; Garcia et al. 1994). In this paper, we 
present results on decomposing flare components into each of 
the two phases in the VUV and the TSI for the five flares 
observed by TIM. We also use a VUV empirical model as a 
proxy for undetected solar flare TSI variations. Our updated 
algorithm is used to revise values in the study by Emslie et al. 
(2012). In Section 2, we discuss the data and empirical model 
used for this study. In Section 3, we describe the time series 
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Table 1 

TIM Observed Flares (in 10 30 erg) 


Date 

Time 

Class 

Location 

TSI-Tot 

TSI-Imp 

TSI-Grad 

VUV-Tot 

VUV-Imp 

VUV-Grad 

2003 Oct 28 

9:51 

X17 

E08 S16 

626 

87 

539 

112 

7.0 

106 

2003 Oct 29 

20:37 

X10 

W10S17 

346 

44 

302 

49 

5.2 

43 

2003 Nov 4 

19:29 

X28 

W83 S19 a 

290 

32 

258 

89 

2.8 

86 

2005 Sep 7 

17.17 

X17 

E77 Sll a 

436 

6.2 

430 

76 

2.6 

74 

2006 Dec 6 

18:29 

X6.5 

E63 S06 

76 

53 

23 

20 

2.1 

18 


Notes. The five X-class flares with detected TIM TSI signatures from solar cycle 23. This list includes total bolometric energy 
release derived from the TSI time series fits. Also included are the VUV (1-190 nm) energy estimates obtained by integrating 
the FISM flare profiles. Both energy sets are decomposed into their impulsive and gradual phases. 
a Designates a limb flare (flare with an east-west position larger than 70). 


fitting procedure used to fit the empirical model to the observed 
TSI data. In Section 4, we explain how these results are used in 
estimating the total X-class solar flare radiative output during 
TIM eclipse periods. In Sections 5 and 6, we discuss results and 
future work, respectively. 


2. DATA 

2.1. Total Irradiance Monitor 

Total solar irradiance data was obtained from NASA’s SOlar 
Radiation and Climate Experiment ( SORCE ) mission by the 
TIM. TIM has been taking solar data since SORCE was launched 
in 2003 January (Kopp et al. 2005a). This instrument sets the 
two criteria in determining which flares are to be studied for 
this analysis: (1) flares that have occurred since 2003 January 
and (2) flares that have an observed signal response in the TSI. 
Table 1 lists the five flares that met these criteria. The SORCE 
spacecraft is in low Earth orbit, which takes about 90 minutes 
to complete. This orbit leaves TIM with roughly 45 minutes 
windows of solar observing when not inside Earth’s shadow 
(non-eclipse periods). The TIM has a time cadence of 50 s. The 
SORCE TIM is the first instrument to have the ability to observe 
solar flares in the TSI (Kopp et al. 2005b; Kopp & Lean 201 1), 
partly due to this time cadence. The TIM rigorous calibration 
and metrology of all its individual components are designed 
to attain a 100 ppm (0.010%) combined standard uncertainty 
in the TSI (Kopp et al. 2005a). This absolute precision gives 
TIM the sensitivity necessary to observe solar flares. Each TIM 
data point has an uncertainty of roughly otsi = 2.72 mW m 2 , 
corresponding to 2 ppm (1 ppm = 1 360/ 10 6 W m” 2 ). 

2.2. Flare Irradiance Spectral Model 

The Flare Irradiance Spectral Model (FISM; Chamberlin 
et al. 2007, 2008) was used to determine the solar irradiance 
in the VUV wavelengths from 0.1 to 190 nm. FISM is an 
empirical model that estimates the solar irradiance in the VUV 
wavelengths with 60 s temporal resolution and 1 nm spectral 
resolution from 0.1 to 190 nm. These qualities allow the 
inclusion of variations due to solar flares, solar rotation, and solar 
cycle variations. FISM utilizes data from the Thermosphere 
Ionosphere Mesosphere Energetics and Dynamics (TIMED) 
Solar Extreme-Ultraviolet Experiment (SEE; Woods et al. 
2005) and the Upper Atmospheric Research Satellite ( UARS ) 
SOlar STellar Irradiance Composition Experiment ( SOLSTICE ; 
Rottman et al. 1993). FISM also includes the derived IP and GP 
for solar flares by using the irradiance and the time derivative 
from the GOES XRS 0. 1-0.8 nm channel. We analyzed 13 solar 
flares of solar cycle 23 in the VUV consistent with the Neupert 


effect (Neupert 1968; Equation (1)). The version of the Neupert 
effect that we employ takes only the positive values of the soft 
X-ray flux time derivative, F SXI , which agrees with the hard 
X-ray flux, E\ lxr , 


d(F s „) 

dt 


Fhxr- 


(1) 


This relation holds when the magnetic loop plasma is heated by 
the influx of accelerated particles and cooling by radiation and 
conduction are negligible (Benz 2008). 


3. DATA ANALYSIS 

The total radiative energy in the VUV wavelengths of each 
flare observed by GOES XRS can be modeled with FISM. The 
FISM formulation calculates a daily flux component separately, 
along with the individual flux contributions from the IP and GP. 
Therefore, contributions from both the IP and GP are already 
known and separate in FISM. The fluence (joule m -2 ) for each 
phase is computed by integrating the irradiance (W m 2 ) of the 
IP and GP over their respective temporal durations of the flare 
and over the VUV wavelengths (0. 1-190 nm) as in Equation (2): 

P tend P '‘'-max 

Evuv = / / / Jfism, ip+gp(^-i t) dh dt , (2) 

7 'start --mi r 

where f start and f cn(1 are the start and end times of each flare phase, 
respectively. The low and high wavelength bounds are 7. mm and 
k max , respectively. The irradiance contributions computed here 
are solely from the flare enhancement and not from the rest of 
the Sun. After the fluence is obtained. Equation (3) gives the 
conversion,/, in units of erg m 2 joule” 1 , which must be applied 
to convert from irradiance units of W m 2 at the Earth to erg at 
the Sun by multiplying 

/ = (1 AU) 2 2tt x 10 7 

/= 1.406 x 10 30 , (3) 

where AU is the astronomical unit in meters. This conversion 
assumes a uniform angular distribution of the flare energy 
release. A similar calculation was performed in Woods et al. 
(2006). This algorithm gives an approximation of the energy 
release in the VUV wavelengths during each phase of the flare 
with similar accuracies to those of FISM flare estimates. The 
VUV radiative energies for all flares in this paper are given in 
Tables 1 and 2. 

An amplitude fitting routine employing the FISM model was 
developed to fit the time profile of flares observed in the TIM 
TSI to obtain the bolometric radiative energy output due to 
solar flares. The TIM TSI data was fit to the scaled IP and GP 
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Figure 1. TSI background subtraction for the 2003 October 28 flare and the 2006 December 6 flare that was not in Woods et al. (2006). The solid line is the linear fit to 
the fluctuating background. The diamonds are the TIM TSI data with error bars from the instrument uncertainties only. The dashed line is the GOES XRS 0. 1-0.8 nm, 
1 minute time cadence data. 


Table 2 

TIM Reconstructed Flares (in 10 30 erg) 


Date 

Time 

Class 

Location 

TSI-Tot 

TSI-Imp 

TSI-Grad 

VUV-Tot 

VUV-Imp 

VUV-Grad 

2002 Apr 21 

0:43 

X1.5 

W84 S14 a 

112 

3.8 

108 

26 

1.0 

25 

2002 Jul 23 

0:18 

X4.8 

E72S13 3 

102 

7.2 

95 

24 

1.9 

22 

2003 Nov 2 

17:03 

X8.3 

E56N14 

148 

43 

105 

47 

4.2 

43 

2005 Jan 15 

0:22 

X1.2 

E14N8 

8.4 

2.1 

6.3 

2.8 

0.2 

2.6 

2005 Jan 19 

8:03 

XI. 3 

W5N15 

64 

16 

49 

22 

1.5 

20 

2005 Jan 20 

6:36 

X7.1 

W61 N14 

218 

40 

178 

77 

3.9 

73 

2006 Dec 5 

10:18 

X 9.0 

E79 S07 a 

166 

19 

147 

39 

5.0 

34 

2006 Dec 13 

2:14 

X3.0 

W23 S05 

113 

37 

76 

35 

3.6 

31 


Notes. The eight X-class flares with reconstructed TSI values based on the FISM model flare profiles. The reconstructed 
TSI impulsive and gradual phases energies are obtained by scaling the VUV impulsive and gradual phase energies as 
described in Section 4. 

a Designates a limb flare (flare with an east-west position larger than 70). 


FISM VUV curves to acquire the energy release in each phase. 
The sum of the energies from each phase gives the total energy 
from the flare in the TIM TSI. The GOES XRS 0. 1-0.8 nm 
channel was used as a temporal proxy for the GP, and the 
time derivative of this channel is used as the temporal proxy 
for the IP. This procedure is the application of the Neupert 
effect. The TIM TSI flare signal is isolated by subtracting the 
background TIM TSI irradiance. A linear relationship is found 
from the irradiance values before and after the flare profile, and 
the resulting background fit is subtracted off (see Figure 1). 
To analyze the TIM TSI flare signal, the GOES XRS (from 
the FISM VUV), IP and GP proxies are scaled such that their 
sum matches the time signature of the background subtracted 
TIM TSI data; see Equation (4), Figure 2, and Figure 3. This 
fit is accomplished by two separate algorithms, an analytical 
and a numerical one. The subscript, E in the definitions refers 
to time in the following expression. The analytical algorithm 
solves for the scaling coefficients of the IP /,. with coefficient 
a, and GP G f , with coefficient /;, of the time profiles exactly by 
a least squares calculation (Bevington & Robinson 2003; Press 
et al. 2007) so that the difference between the observed data and 


model function, Z f , is minimized: 

Z,=aI,+bG t . (4) 

The values obtained by mathematical definition will be the best 
coefficients if the errors in the data are Gaussian distributed 
and the errors are relatively small compared to the data value. 
Our error bars (solely from instrument precision and background 
subtraction), as displayed in Figures 2 and 3, are small compared 
to their corresponding values. The data value errors are not 
necessarily Gaussian distributed, but this merit function still 
gives a good estimation of the best scale factors. 

The values obtained by this methodology can be positive or 
negative, where negative scale factors for irradiance are non- 
physical. A negative scale factor was encountered for the IP 
of the flare on 2005 September 7, due to a very small IP. In 
addition to the analytical algorithm, a numerical minimization 
of chi-square was conducted by the exploration of parameter 
space along the minimum “bowl” for a predefined grid of values 
for the a and b coefficients. This method is limited by the choice 
of gridded values for a and b that are explored. Generally, the 
analytic and numerical scale factors are consistent with each 
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Figure 2. Background subtracted data for four of the five X-class flares: X17 flare on 2003 October 28, X10 flare on 2003 October 29, X28 flare on 2003 November 4, 
and the X17 flare on 2005 September 7. The diamonds are the background subtracted TSI data with error bars comprised of instrument uncertainties added in quadrature 
with the linear background fit uncertainties. The dashed line is the scaled FISM gradual phase profile. The dashed-dot-dot-dotted line is the scaled FISM impulsive 
phase profile. The solid line is the sum of the scaled impulsive and gradual phases. The vertical dotted lines are the timeframes of which the TSI data was fit to the 
FISM time profiles. 



Time (Hours) 

Figure 3. Background subtracted data for the X6.5 flare on 2006 December 6 not 
included in Woods et al. (2006). The diamonds are the background subtracted 
TSI data with error bars comprised of instrument uncertainties added in 
quadrature with the linear background fit uncertainties. The dashed line is the 
scaled FISM gradual phase profile. The dashed-dot-dot-dotted line is the scaled 
FISM impulsive phase profile. The solid line is the sum of the scaled impulsive 
and gradual phases. The vertical dotted lines are the timeframes of which the 
TSI data was fit to the FISM time profiles. 


other. The analytic algorithm for scale factors were used for 
all flares in this paper, except for the 2005 September 7 flare 
for the aforementioned reason. The scaled GOES time series 
were integrated with respect to time, and the conversion factor 
(Equation (3)) is applied to obtain the TSI contribution for each 
flare phase at the Sun. The radiative energies for the five flares 
with an observable TIM TSI signature are given in Table 1 . 

4. RECONSTRUCTED TSI FLARE ENERGY RELEASE 

The radiative output from X-class solar flares have large 
contributions from the VUV wavelength emissions, ~41% 
for the GP of center flares and ~23% for the GP for limb 


flares. We take limb flares to have an east-west location >70°. 
Using the five flares seen in the TIM TSI and estimated 
VUV emission simultaneously obtained from TIMED SEE 
and GOES, respectively, the TSI flare irradiances during TIM 
eclipse periods can be reconstructed. This method can be used 
to estimate the position-dependent total radiative output from 
flares, which is described in this section. 

The average position-dependent fractional energy contribu- 
tion to the flare from the VUV to the total radiative output is 
found for each flare. The inverse of these values are used to 
scale the VUV energy output for flares not observed due to TIM 
eclipse periods (when TIM is behind the Earth with respect to 
the Sun) to obtain an estimation of the total radiative output. 
This model is still limited by low statistics (five solar flares), 
but as more large flares occur, more robust algorithms can be 
developed. The IP Amp and the GP Be, mil scale factors for the 
time-integrated VUV irradiance used to approximate the time- 
integrated TSI flare output are given in Equation (5) for the IP 
and Equation (6) for the GP. The total number of limb or center 
flares is given by n, the index for each individual flare in each 
category is i, the VUV IP energy is VI, the TSI IP energy is 
T I, the VUV GP energy is V G, and the TSI GP energy is TG: 



(5) 

(6) 


These calculations are done separately for center and limb 
flares. The scale factors for center flares are Amp = 10.3 and 
^Grad = 2.4 for the IP and GP, respectively. Limb flare scale 
factors obtained were Ai mp = 3.8 and Z?Grad = 4.3. 
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The TSI energy estimations are the same order of magni- 
tude as previously published values for X-class solar flares 
(10 31 — 10 32 erg) and have accuracies consistent to the order given 
by the uncertainties. TSI-radiated energies using the previously 
described algorithm vary somewhat from the values given in the 
TIM analysis by Woods et al. (2006) by up to a factor of 1.5, 
where an additional flare on 2006 December 6 is included here 
that occurred after publication of the aforementioned reference. 
This is believed to be due to a more robust fitting routine for both 
the IP and GP scaling, but more importantly, to the improved 
background fitting routine that can easily give a factor of two in 
the final results due to the long tail in the decay phase that flares 
have. The latter is most apparent in the 2005 September 7 flare, 
which is dominated by GP emission. 

Our results are also consistent within the 1 a uncertain- 
ties of the TSI flare contributions derived in the analysis of 
Kretzschmar et al. (2010) using the PMO and DIRAD radiome- 
ters from the VIRGO experiment on the SOlar and Heliospheric 
Observatory mission (Frohlich et al. 1997). As anticipated, the 
GP contains the bulk of the energy released in solar flares as 
compared to the IP (an average of 7.1 times larger), due to the 
GP lasting much longer. Tables 1 and 2 demonstrate variations 
in the radiative output of the flares observed at Earth due to 
center-to-limb variations of optically thick emissions (Worden 
et al. 2001; Chamberlin et al. 2008). Limb flares have an ex- 
pected lower energy output because the XUV wavelengths from 
optically thin coronal emissions dominate the energy release and 
are unaffected. The optically thick emissions are reduced due to 
the extended path length through the solar atmosphere. Flares 
near disk center are more intense than the limb flares in the TSI. 
This center-to-limb correction can also be applied in reverse 
to get the TSI radiated output for a limb flare as if it actually 
occurred on the center of the solar disk. The modeled values are 
consistent with values given for two flares that occurred before 
2003 that were analyzed in Emslie et al. (2004). 

This model can be extended for use on flares lower than 
X-class, which was done in an earlier version of this analysis 
in the investigation by Emslie et al. (2012). In that study, the 
reconstructed TSI was computed for the majority of X-class 
flares of solar cycle 23 that were not observed by TIM due 
to eclipse periods or flares that occurred before the launch of 
TIM in 2003. This matured version provides a more accurate 
estimate of the radiated energy for the energetic comparisons of 
SEEs. The revised values given here are, in general, larger by 
roughly a factor of two than those from the Emslie et al. (2012) 
project due to the modified processes to subtract the highly 
variable background and, equally as important, the time interval 
in which the fit is performed and the selection of adequate start 
and end times for temporal integration of the TSI data. 

These newer values are consistent with the conclusions 
derived in the comprehensive study of Emslie et al. (2012). The 
bolometric radiated energy is still over an order of magnitude 
less than the stored non-potential (free) magnetic energy in 
the studied active regions. This still agrees with the belief that the 
available magnetic energy is more than sufficient to power the 
flare components. The bolometric radiated energy is still less 
than the coronal mass ejection kinetic energy. Moreover, the 
bolometric radiated energy is comparable to the lower limits 
of the flare-accelerated particles, where the latter could be 
up to an order of magnitude greater due to the lower energy 
cutoff being masked by the thermal plasma radiating in soft 
X-rays. The conclusion persists that there is sufficient energy 


in accelerated particles to power the radiative energy release. 
Finally, the correlation between the soft X-ray plasma and 
bolometric energy is preserved. 

6. DISCUSSION AND FUTURE WORK 

A full in-depth knowledge of the energy distribution of solar 
flares across the entire electromagnetic spectrum is difficult 
to obtain due to lack of simultaneous observations by the 
instruments at the time of analysis. These modeled results 
are the best estimations currently achievable, but they are 
limited to the current TIM TSI data and spectrally to the VUV 
region. We have broken down the radiative contribution of the 
VUV (0.1-190 nm) wavelengths to the total radiative energy 
release for the largest solar flares, partitioned them into two 
separate phases, and have implemented an empirical model for 
bolometric flare energy reconstruction. In the future, more steps 
will be taken to pursue increased knowledge of the spectral 
energy distribution of flares along with IP and GP composition. 

The authors thank the Laboratory for Atmospheric and Space 
Physics (LASP) for their funding and support, as well as support 
from the NSF REU program. 

APPENDIX 

For completeness, our analytic fitting method is described 
here. The subscript, t, in the definitions refer to time in the 
following expressions. The analytical algorithm solves for the 
scaling coefficients of the IP /,, with coefficient a, and GP 
G t , with coefficient b, of the time profiles exactly by a least 
squares calculation similar to the one in Press et al. (2007). The 
algorithm is as follows. For a model function, Z,. that is linear 
in the parameters a and b (Equation (Al)): 


Z,=aI,+bG,. (Al) 

The coefficients that give the optimal approximation to the 
actual observed TSI data, D t , with N discrete data points, and 
propagated uncertainties <j t , can be found by composing a 
chi-square merit function, differentiating this function with 
respect to parameters a and b separately, setting them equal 
to zero and then solving these two coupled equations for the 
coefficients a and b (Equations (A2)-(A1 1)): 


2 (<z, b) = Y] 


i=i 


D t — al, — bG t 


(A2) 


Partial differentiate with respect to a and b, set equal to zero for 
two separate equations: 


d(x 2 (a, b)) 

da 



(A3) 


9(x~(a, b)) 
db 


= E-2 


i=i 


D, — al, — bG t 


G, = 0. (A4) 


We use the substitutions, 



(A5) 
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s - V 0,1 

bDl ~ 

i= i a > 


(A6) 


S/c = X! 


/,G, 


t= 1 Gt 


(A7) 


5 dg = X! 


D,G, 


t= 1 


(A8) 


and then solve for a and i>. 

_ SdiSgg ~ SdgSig 
SiiSgg ~ ( Sig ) 2 

SiiSdg — SdiSjg 
SiiSgg — (Sig) 


b = 


,2 ' 


(A9) 


(A10) 


(All) 


This gives the impulsive and GP scale factors, a and b, 
respectively. 
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